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We investigate the propagation of microwave photons in a one-dimensional open waveguide interacting with 
a number of artificial atoms (qubits). Within the formalism of projection operators and non-Hermitian Hamil¬ 
tonian approach we develop a one-photon approximation scheme for the calculation of the transmission and re¬ 
flection factors of the microwave signal in a waveguide which contains an arbitrary number N of non-interacting 
qubits. We considered in detail the resonances and photon mediated entanglement for two and three qubits 
in a chain. We showed that in non Markovian case the resonance widths, which define the decay rates of the 
entangled state, can be much smaller than the decay width of individual qubit. It is also shown that for identical 
qubits in the long wavelength limit a coherent superradiant state is formed with the width being equal to the sum 
of the widths of spontaneous transitions of N individual qubits. The results obtained in the paper are of general 
nature and can be applied to any type of qubits. The specific properties of the qubit are only encoded in the two 
parameters: the qubit energy Q and the rate of spontaneous emission T. 
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I. INTRODUCTION 


One-dimensional (ID) waveguide-quantum electrodynam¬ 
ics (QED) systems are emerging as promising candidates for 
quantum information processing motivated by tremendous ex¬ 
perimental progress in a wide variety of solid state systems 
with imbedded artificially designed atoms- qubits I l! ;T'|. Con¬ 
fining the microwave field in reduced dimensions such as ID 
waveguide and taking account of the enormous dipole mo¬ 
ment of artificial atom the photon qubit interaction can be 
strongly enhanced as compared with open 3D space HU). 
In recent years one of the basic type of these physical systems 
has been realized in solid state setups where qubits were on 
chip coupled to microwave cavities H. An important advan¬ 
tage of these systems is that the qubits can be placed within 
the photon field confined in a microwave cavity at fixed pre¬ 
determined positions at separations on the order of relevant 
wavelength. Moreover, unlike the real atoms, qubits are in¬ 
trinsically not identical due to technological scattering of their 
parameters. It is also important that the excitation energy of 
every qubit in a chain can easily be adjusted by external cir¬ 
cuit. 

The experimental investigation of these systems is based 
on the measurements of the transmitted and reflected signals 
with their properties being dependent on the quantum states of 
every qubit in a waveguide. Up to now there are known only 
several experiments_with a single superconducting qubits in 
ID open space |!8- IS] and one experiment with two transmon¬ 
type qubits in a waveguide. 0. 

For solid-state quantum information processing it is inter¬ 
esting to study ID waveguide systems having more than just 
one qubit. A key point here is whether the multi-qubit system 
could display a long lived entanglement necessary for imple- 
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mentation of quantum algorithms. The entanglement is also 
necessary for quantum error correction which requires at least 
three qubits in the chainfl7j. 

Such multi-qubit systems exhibit, in general, non Marko¬ 
vian behavior: the interaction between qubits is not instanta¬ 
neous, hence, the retardation effects have to be included. The 
manifestation of these effects is that the resonances (energies 
and their widths) of the qubit system become dependent on 
the frequency of incident photon. In this case a master equa¬ 
tion for the density matrix p of the qubits cannot be written 
in Lindblad form. A Markovian approximation corresponds 
to long wavelength limit, kd « 1, where k is photon wave 
vector, d is a distance between neighbor qubits. In this case we 
may neglect the retardation effects and assume that the qubits 
interact instantaneously. 

Recent experiments with superconducting qubits showed 
that the photon mediated interaction between distant qubits 
can lead to the creation of two-qubit lfl6L H~8h and multi-qubit 
entanglement JT^-jll]] . 

Theoretical calculations of microwaves transmission in ID 
open waveguide with a qubit placed inside are being per¬ 
formed in a configuration space ll22l - l26ll or by the input- out¬ 


put formalism 1127112811 . These methods are physically sound 
but they become very cumbersome if we try to find solutions 
for several or more qubits in a waveguide. While the trans¬ 
mission for a single two level atom in ID open waveguide has 
long been known ||22i l23ll . the analytical expressions for the 
transmissions for two qubits and for symmetrical arrangement 
of three identical qubits have been published quite recently 


For N identical equally spaced qubits the transmission can 
be found analytically with the help of the method borrowed 
from the physics of crystals with translational symmetry |29i]. 
However, in general, for N non identical randomly spaced 
qubits a simple analytical procedure does not exist. 

In the present paper we propose a matrix formalism for 
the study of a one- photon transport in ID open waveguide 
filled, in general, with a number of not identical and arbitrary 
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spaced artificial atoms. Similar idea has been suggested for 
the study of photon transport in the coupled resonator opti¬ 
cal waveguides 0. Our approach is based on the projec¬ 
tion operators formalism and the method of the effective non- 
Hermitian Hamiltonian which are the powerful tools to deal 
with a Lippman-Shwinger scattering problem. It is different 
from usual resolvent method of solving Lippmann-Schwinger 
equation which was demonstrated for two qubits in (3j|. 

The method we use here has originally been developed for 
the description of nuclear reactions |0 ;3j] with many later 
applications for different open mesoscopic systems ranging 
from universal conductance fluctuations D4|1 to electron trans¬ 
port through ID solid state nanostructures [3H[36|](see review 
paper (37;] and references therein). 

For general N qubit case our technique allows us to easily 
include the cases of non identical qubits and/or with unequal 
spacing when the translation symmetry is absent. It is very 
important for artificial atoms with inevitable technological 
spreading of parameters with the possible individual tuning 
of qubit resonance energies. Additionally, a direct exchange 
interaction between nearest neighbor qubits can be easily in¬ 
corporated into the scheme of this technique. The influence 
of this interaction between two and N superconducting flux 
qubits on the photon transmission and entanglement has been 
studied in the papers H^-0]. 

With the aid of our technique we study in detail the one- 
photon microwave transport for one, two and three qubits 
imbedded in a waveguide. We considered in detail the res¬ 
onances and photon mediated entanglement for two and three 
qubits in a chain. We show that for N identical qubits in 
the long-wavelength limit a coherent superradiance state is 
formed with the width being equal to the sum of the widths 
of spontaneous transitions of N individual qubits. 

The paper is organized as follows. In the Section II we de¬ 
fine the model Hamiltonian of N noninteracting qubits imbed¬ 
ded in a microwave resonator. In Section III we describe 
in detail a projection formalism and effective non-Hermitian 
Hamiltonian approach in application to the photon transport 
in ID waveguide. The application of the model Hamiltonian 
to the derivation of the general expressions for the transmis¬ 
sion and reflection coefficients for N qubits in a waveguide 
is performed in the Section IV. The Section V is devoted to 
a detailed investigation of the microwave transport for one, 
two, and three qubits in a waveguide. In this section we give 
not only the analytical expressions for the transmission and 
reflection factors for two and three qubits in general case, but 
we investigate in detail the energy spectrum of these systems 
and their resonances in non Markovian case, which is auto¬ 
matically included in our theory, since the quantity kd explic¬ 
itly enters the analytical expressions. We also study a photon 
mediated entanglement for two and three qubit systems. The 
results of this section are important for three qubits experi¬ 
ments, since to our knowledge there are no ID open space 
experiments with three qubits in a waveguide. In the conclu¬ 
sion to this section we briefly analyze the general case of N 
qubits. 


II. THE MODEL HAMILTONIAN 

We consider a microwave ID waveguide resonator with N 
qubits imbedded at the fixed positions Xi. The Hamiltonian of 
the system reads: 

H = H ph + H qb + H int (1) 

where 

H ph = a k (2) 

k 

is the Hamiltonian of photon field, 

H 'd>= e:, h i<> (3) 

is the Hamiltonian of N noninteracting qubits, where 

H\ b = (4) 

is the Hamiltonian of the individual i-th qubit with the excita¬ 
tion frequency O,. 

The interaction of the qubit chain with the photon field is 
given by Hamiltonian: 

Hint = EL M4 e ~ ikXi + a k e ikXi )a^ (5) 

k 

where A, is the qubit-photon interaction strength, are the 
qubit positions relative to the waveguide center, xq = 0. 

III. PROJECTION FORMALISM AND EFFECTIVE 
NON-HERMITIAN HAMILTONIAN 

1. Projection operators formalism from the formal point of view 

As the projection operators formalism and effective non- 
Hermitian Hamiltonian approach are not common in the field 
of quantum optics, here we briefly describe the essence of this 
method omitting its rigorous justification which can be found 
in the corresponding literature (see review paper 0 and ref¬ 
erences therein). 

It is always possible to formally subdivide the Hilbert space 
of a quantum system with the Hermitian Hamiltonian H into 
two arbitrarily selected orthogonal projectors, P and Q, which 
satisfy the properties of completeness 

1 = P + Q (6) 

and orthogonality 

PQ = QP = 0 (7) 

From (|6| and £7} it also follows: 


PP = P, QQ = Q 


( 8 ) 
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With the help of the completeness © we can divide the 
solution of the stationary Schrodinger equation, 

H = E® (9) 

in two parts, 

^ = + = + ( 10 ) 

and rewrite ((9} in the following form: 

(P + Q)7T(P + Q)('Tp + 'T q ) = P('Tp + 'T q ) (11) 

Since, in virtue of (|7}. P4 /q=0, Q y \> p = 0 we rewrite (fTTb as 
follows 

(Hpp + Hqp)^ p + (Hqq+Hpq)^q = E('&p + '&q) 

( 12 ) 

where 


the discrete states. We also assume that Hamiltonian Hpp 
is diagonal in subspace P. In order to avoid the singularities 
emerging when Hpp has eigenvalues at real energy E, it has 
to be considered as a limiting value from the upper half of 
the complex energy plane, E + = E + ie. With this rule, 
those states of subspace Q which will turn out to be coupled 
to the states in subspace P will acquire the outgoing waves 
and become unstable. Then, for this scattering problem the 
effective Hamiltonian 06} becomes non Hermitian and has to 
be written as follows: 

Heff(E) = Hqq + H QP E _ H 1 pp+i£ H PQ (17) 

In this case the equation 03} defines the resonance energies of 
the Q- system which lie in the low half of the complex energy 
plane, E = E — iHT and are given by the roots of the equation 

det (E — H e ff) = 0 (18) 


H qq = QHQ , Hpp = PHP, Hqp = QHP, H PQ = PHQ 

The equation 03} is equivalent to the Schrodinger equation 

©. 

Multiplying 03} from the left by projectors P and Q we 
obtain two coupled equations for 4/ p and Tq: 

{Hpp - E)^ P = -Hpq^q (13) 

{Hqq - E)^q = -Hqp'Hp (14) 

If we eliminate from ( IT3l > or 04} one subspace of states, we 
can obtain an equation for a part of the wave function P0} . 

For example, if we eliminate from 03} the P-subspace, the 
equation for the wave function in Q- subspace takes the form 

H eff (E)* Q =E^q (15) 

where the energy dependent effective Hamiltonian 

H eff{E ) = Hqq + Hq P ^ H P Q (16) 


The imaginary part T of the resonances describes the decay of 
Q- states due to their interaction with P- states. 

In the framework of projection formalism we can find from 
03}, 03} the wavefunction of the whole system 4', a solution 
of Shrodinger equation (|9}, in terms of the operator which acts 
on the initial state, | in), which contains continuum variables 
and satisfies the equation Hpp\in) = E\in), where E is the 
same as in (|9}. Then, the formal solution of 03} can be ex¬ 
pressed in the following form 

H/p = \in) + — -—-—- Hpq^q (19) 

E — Hpp + ie 

Substituting this expression in r.h.s. of 04} we obtain 

= f \ R Q p \ in ) (2 °) 

E - Heff 

where fT e // is given by its non Hermitian form 07} . As a 
final step, we substitute Tq from ( 120} into r.h.s. of 09} and 
combine these two equations to obtain the expression for the 
state vector of the Shrodinger wavefunction T [41] 


projects Hilbert space on the Q subspace. 

The second term in 06} describes multiple excursions to 
the class P with return to the class Q. 

It should be noted that while the energy E in 03} is the 
same as in Schrodinger equation ([9}, the equation 07} is not 
equivalent to (0: the effective Hamiltonian 06} is energy de¬ 
pendent , so that the eigenvalue E enters this equation in a 
complex way, and wavefunction Tq is not egenfunction for 
E, however it can be written as a linear superposition of the 
state vectors from Q- subspace. 


|'T) = | in) + 


1 


1 


E — Hpp + ie 


-H 


E ~ H eff 
1 


PQ 


E-H, 


eff 


Hqp | in) 
Hqp | in) 


( 21 ) 


In fact, this expression is nothing more than a decomposition 
PO} . The last term in ([31} is the part of 4/p, which describes 
to all orders of Hqp the evolution of initial state |in) under 
the interaction between P and Q subspaces. 


2. Application to the scattering problem 

Keeping in mind the scattering problem we assume that Q 
subspace consists of discrete states, and P subspace consists 
of the states from continuum, however, it may also contain 


3. One photon scattering 

In the one photon approximation there are two possibil¬ 
ities: either one photon is in the waveguide in the state 
|lfc) = | k) and all qubits are in their ground states \gi) with 
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a corresponding state vector \gi,g2 . gN,k), or no pho¬ 

tons in the waveguide, | 0 fc) = | 0 ), with ?-th qubit being 
excited and N — 1 qubits being in their ground states. In 
this case the system is described by N vectors of the type 

|si> ■■gi-i,e i ,g l+ i,...gN , 0 ). 

In order to simplify the notations we will use throughout 
the paper the following concise forms for state vectors: 

\k) = \gi ,92 . gN,k) (22) 


where Rm,n is the matrix inverse of the matrix ( m(E — 

H ef /)\n): 

Rm,n = {m\— -— 77 —I n) (29) 

- tieff 

The second term in (l28l > is the wavefunction Tg of a qubit 
system modified by its interaction with a photon field: 

N 

’Fq = ( to | Hqp \k) (30) 

n,m=l 


|^) — ”9n— 1 ? Qn+h "-9Ni 0) (23) 


with the orthogonality relations 

( n\m) = S nm 

(n\k) = 0 

(k\k') = -5(k-k'), 

J-j 

where L is the waveguide length. 

In these notations the initial state is just the state d 
(| in) = |fc}) with the energy 

h N 

E = E k = hio k - - ^2 (24) 

Z i=l 

where ujk is the frequency of incident photon. 

Hence we take the projection operators as follows: 

+oo 

P = J2 k I k) (k\ = A J dk | k) (k\ (25) 


N 

Q = l n ) ^ n l (26) 

n=1 

Then the matrix elements of the effective Hamiltonian (fTTI i 
in subspace Q can be written as 


In more general context the expression (l30l > describes the en¬ 
tanglement between qubits due to their interaction with a pho¬ 
ton field. 

From (l30l > we can also find the probability for the n- th 
qubit to be in excited state: 

N 

(n\^ Q )=J2Rn,m(m\H QP \k) (31) 

m =1 

The photon wavefunction in configuration space is ob¬ 
tained by multiplying (l28l > from the left by the vector (a;| = 
(x,g N ,gN-i, . gi\: 

^n(x) = (a.jT') = e ikx + 

£ E Jdq Ek l'E q+ie {q\ H P q 1 n) R n ,m <m| H QP |fc) 

n,m—1 

( 32 ) 

where we have used the definitions ( x\k) = e lkx and ( x\n) = 
0 . 

The wavefunction (|32] > is a superposition of the incident 
wave and the wave which results from the virtual transitions 
between qubits and photon field in the resonator. We will see 
below that this superposition leads to the destructive interfer¬ 
ence when the frequency of incident photon is equal to the 
excitation frequency of any qubit. In this case the transmitted 
signal outside the qubit array is equal to zero. 

IV. THE APPLICATION OF PROJECTION FORMALISM 
TO THE MODEL HAMILTONIAN 


( m\H eff \n) = (m\H\n) + 


+oo 

L_ [ d (m\H QP \q ) (g| H PQ \n) 
2n J q E k - E q + ie 


(27) 


— OO 


In the basis of Q- subspace vectors the full wavefunction 
d can be written as 

N 


W = |fc) + E \n)Rn,m (m\ H qp \k) + 

n,ra=l 

h E I d( l E k -E q +ie (gl Hp Q I 71 ) Rn ,m {m\ Hq P \k) 

n,m =1 


Here we apply the model Hamiltonian ([TJ from the Section 
mho the general expressions found in Section|III] First we cal¬ 
culate the matrix elements of the effective Hamiltonian (l27l >. 
The first term in rhs of (l27l> reads: 

(m\ H | n) = £ m 5 m ,n (33) 

where 

£m = 

It is also not difficult to calculate the matrix elements in rhs of 
equation (1271 : 




N 

E 




(34) 


(28) 


(m\HQ P \k) = X m exp(ikxm) 


(35) 
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Then, the the second term in rhs of (IT 1\ can be written as 

(36) 

where 


J(a 


+ OO 

i ? *^n) — ^ & 


dq 


exp (iq{x m - x n )) 
Ek — Eq + is 


(37) 


It is shown in the Appendix that 


where the transmission and reflection coefficients are as fol¬ 
lows 


N 

t N = l-ih J2 (r m r n ) 1/2 e ikXm R m , n e~ ikx " (45) 

m,n =1 


N 

r N = -in Y (TmT n ) 1/2 e ikx ™Rm,ne ikXn (46) 

m,n= 1 


J( a 


2in 

hv g 


e ik\d mn \ 


(38) 


The conservation of the energy flux requires the additional 
condition for t and r : 


where d m n = x m — x n , and k is related to the physical fre¬ 
quency u> of incident photon, k = uj/v g , where v g is the group 
velocity of the photon wave in a waveguide. 

Finally, the effective Hamiltonian (IT7b can be written as fol¬ 
lows 


\t N \ 2 + \r N \ 2 = l (47) 

The expressions (l45l > and (l46l> are of general nature and they 
form the basis for the calculation of microwave transmission 
and reflection in particular cases. 


(m\H eff \n) = £ m 5m,n — ih{T m T n ) 1 l 2 e ik ' [dmrl ' { (39) 

V. TRANSMISSION, REFLECTION, AND PHOTON 

where we define the halfwidth of spontaneous emission MEDIATED INTERACTIONS IN THE QUBIT SYSTEM 


T A 2 

r m = -r^ (40) 

n 2 v g 

Throughout the paper we will use T for the halfwidth of 
resonance line. 

The photon wavefunction (l32l > for our model follows from 
(ED and ®: 


A. One qubit in a waveguide 

In this case, in subspace Q there is the only vector |1). We 
also assume that the qubit is located at the point x = 0. From 
(l45l > and (l46l > we obtain: 


N 

*n{x) = e ikx -ih Y C r m Tn) 1/2 e ikXm R m , n e ik \ x - x 

m,n= 1 

(41) 

where the matrix Ii m .„ is defined in (refRmn). 

Finally, for our model we write down the qubits’ wavefunc¬ 
tion 'hq and the probability for the //-th qubit to be in excited 
state: 


N 

^Q= Y I n)\mRn,mE kXm (42) 

n,m =1 

N 

<n|tf Q > = ^ A m Rn, m e ikx ™ (43) 

m—1 

We assume that all qubits are arranged in the array from left 
to right, so that x\ is the position of the qubit at the left end 
of the array and xn is the qibut’s position at its right end. In 
this case the photon wavefunction (OTT ) outside the array can 
be written as: 


t± — 1 — ihV Rn 


ri = —ihTRu 


where 


-Rn = <1| 


E - II, 


eff 


-ll} = 


E-(l\H eff \l) 


(48) 


The running energy E in (l48l > is the energy of incident pho¬ 
ton plus the energy of the qubit in the ground state, E = 
hu — Md/2. From ( l39t we also have: 

Hence, for t and r we finally obtain: 


LO — 

oj — n + *r 


(49) 


r i 


-*r 

w — f 2 + *r 


(50) 


^n(x) 


tNe lkx (x > xn) 
e ikx +r N e~ ikx {x < xx) 


(44) 


The plots of microwave transmission and reflection are shown 
in Fig|T| At resonance the signal transmission is zero. 
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FIG. 1: Color online. Transmittance T = \t\ 2 and reflectance R = 
|r| 2 as functions of u>/Sl for one qubit in a waveguide. The width of 
resonance is 2T. 


The expressions (|49] , ( I50| i coincide with those obtained in 
ll23l l where one qubit problem has been solved in a configura¬ 
tion space. 

From ( l49b and (l50l > we see that t—r = 1. Unlike the general 
condition (147] it is valid only for one qubit and reflects the 
continuity of the wavefunction (144] at the point x = 0. 

The probability for the qubit to be excited is given by (143] 
with n = 1 : 


{1\*q)=\Rii = ± -= (51) 

nut — il + tL 

B. Two qubits in a waveguide 

1. Spectral properties of effective Hamiltonian 

Here we consider the first nontrivial example that exhibits 
superradiant transition: the two noninteracting qubits in a 
waveguide . The qubits are positioned at the points x\ = 
— d/2 and x 2 = +d/ 2 , respectively, with a distance d be¬ 
tween them. The Q- subspace is formed by two state vectors 
|1) = |ei,<72,0) and |2) = |gi,e 2 ,0). According to (|39] the 
matrix of effective Hamiltonian is as follows: 


Heff = 


e - ihT x -ihs/T^e ikd 
-ihs/YiT^e ikd -e - ihV 2 


(52) 


where e = |(fii — H 2 ), and Tj, (i = 1 , 2 ) are defined in 

@. 

From the matrix (152] we can find the complex energies 
of the Q- system from the equation (fl 8 ] . where E = hui — 
(Hi + fi 2 ). For two qubit case the equation ([!¥] gives two 
poles in the complex u> plane as the function of physical fre¬ 
quency u. 



FIG. 2: The relation between real and imaginary part of the roots in 
the complex ui plane. For every kd there are two roots which lie on 
the circle at the opposite points 


From this expression it follows that the positions of reso¬ 
nances and their widths depend on the frequency u> of inci¬ 
dent photon (k = uj/vg). This is a common feature of non 
Markovian behavior if the number of qubits is more than one. 

For identical noninteracting qubits Hi = fl 2 = S2, Fi = 
r 2 = r we obtain from (153] 

w = n - ir ± iTe ikd (54) 

In the complex uj plane the roots are as follows 

Reu; = E = H T sin kd (55) 

Imuj = T = —I (1 cos kd) (56) 

From ( |55] and ( |56] we obtain the relation between the real 
and imaginary part of the roots 

(£-n ) 2 + (f + r) =r 2 (57) 

It is remarkable that this relation does not depend on the k, 
i. e., on the running frequency uj. In the plane (T, E) all 
solutions of eq. ( 1571 ) lie at the circle centered in the point 
—F, il with radius equal to F. In the long wavelength limit 
(kd << 1) we obtain from (154] two poles 

uj-x- = H -f r —d — i2T 

V 9 

w_ = H-r— 

V 9 

We see that in this approximation one of the states absorbs 
the width of two qubits. With the increase of F two states 
repel each other. This is also holds if kd is integer multiple of 
7r: one state becomes stationary while the width of the other 
state is 2T. However, this case is valid only for particular 
values of the running frequency uj n = nnv g /d, (n = 1 , 2...). 


u> = 


Hi + n 2 
2 


.Tr + r 2 

*—^— ± 


2. Calculation of microwave transmission 



H 2 + i[r 2 -r 1 ]) 2 -r 1 r 2 e 2ifcd 


(53) 


In order to find transmission and reflection factors t and r, 
it is necessary to calculate the matrix R m ,n, ( m,n = 1 , 2 ) 
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which is the inverse of the matrix (E — 7? e //)m,ra, where the 
physical energy E = hut — /i(ft-| + ft 2 )/2, and the elements 
of the matrix (-ff e //)m,n are given in (f52l) . Direct calculations 
yield for R m ,n the following result: 

p _ 1 fui~n 2 + iT 2 -iVT^e ikd \ 

m ’ n HD 2 {uj) \-is/T^T 2 e ikd ui-fl 1 +iT 1 J 

where 


d 2 (u>) = [w - ft 2 + ?r 2 ] [ u - ni + iri]+rir 2 e i2fcd (59) 

Finally, according to prescriptions in (1451) and < 146b we ob- 
tain t and r in terms of running frequency to: 


^2 


_(tu-fli )(q;-ft 2 )_ 

[w — ft 2 + *r 2 ] [w — fii + «r 1 ] + rir 2 e 2 * fcd 



. {e lkd T 1 [cu - ft 2 - zr 2 ] + e~ ikd T 2 [u - Or + H^]} 

[w — ft 2 + *r 2 ] [w — fii + *Fi] + rir 2 e 2ifcd 

(61) 

For identical qubits we obtain from (l60t . (|6H : 


c) d) 

FIG. 3: Color online. The dependence of transmission (black) and 
reflection (red) amplitudes on the frequency of incident photon, uj/Q 
for different values of kod for two identical qubits. 


f 2 


(c 0 - ft ) 2 

D l 2 d (ui) 


(62) 


2 r[(w 

r 2 = -1 - 


ft) cos kd + r sin kd\ 


where 


Df{u) = (w - ft + iT) 2 + r 2 e 2ikd 


(63) 


(64) 


As is seen from these expressions the form of the transmis¬ 
sion and reflection spectra depend on the inter qubit distance 
d. In the long wavelength limit we obtain from d62l ) and (163 b : 


3. Photon mediated entanglement of two qubits 

For two qubits the structure of the function '\’q ( 142b within 
a subspace of qubit states |1) and |2) is a linear superposi¬ 
tion of the two two-qubit states = a |1) + b |2), where, in 
general, a and b depend on the physical frequency ui. For 
two identical qubits we obtain a general expression which 
describes the frequency dependent entanglement of two two- 
qubit states: 

\ xkdj 2 

*<2 = T^TT ([( w - n - <Te«-]|l> 

U 2 

+ (w-fi)|2» (67) 




ui — ft 

ui — ft + i2T 


In the long wavelength limit kd « 1 the maximally entan¬ 
gled superradiant state which corresponds to a coherent sym¬ 
metric superposition is formed: 


r 2 


-* 2 r 

ui — ft + * 2 r 


( 66 ) 


The expressions ( l65t . ( l66l ) are identical to the one qubit case 
(l49l >. (l50t with the only exception. For two identical qubits the 
resonance width is twice the resonance width for one qubit, 
which is clear signature of superradiance transition which cor¬ 
responds to a coherent symmetric superposition ('Fq)s = 
a (|1) + |2)), whete the quantity a is given in the next sub¬ 
section. 

Below we show several plots of transmission and reflection 
amplitudes for different values of kod, where ko = ft/u g . 
The plots are calculated for two identical qubits from (f62l) and 
(l63l >. The points of the full transmission corresponds to ze¬ 
ros of the numerator of the expression < [63l . Hence, the full 
transmission is observed at the points where the reflection is 
exactly equal to zero. 


(*q)s = -p (I 1 ) + |2» (68) 

ui — ft + 2il 

The transmission and reflection in this case are given by the 
expressions (l65l > and (l66l >. The resonance line of superradiant 
state is directly given as the line of reflection factor (l66l i. 

However, for arbitrary values of kd maximally entangled 
states are formed only for particular values of the frequency 
ui. For example, if kd = uid/v g = mr (n = 1, 2,...) we obtain 
from (l67l > the expression 


I'd = 


A i r 


-n + 2ir 


[(— 1 )"| 1 ) + 12 >] 


(69) 


where ui n = mtv g /d. 

For on resonant excitation (ui = ft) and kod mr, 
where ko = ft/t> g we get from (l67b unentangled state T'q = 
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i (A/r) e ife ° d / 2 |l). in this case we observe a full reflection 
with only the first qubit being excited. 


4. Resonances in two- qubit system 

As it follows from the results of subsection IIII 21 the reso¬ 
nances (their energies and widths) in multi-qubit system are 
given by the roots of equation (fl8l) . The widths of these res¬ 
onances define, in general, the decay rates of the qubit wave- 
function T q (l30l) . 

For two qubits these roots, which are labelled below as 
wi,W 2 , are given in (f53l) . The denominator (1 59 1 > can then 
be written as D 2 (oj) = [w — wi(w)] [w — w 2 ( 0 ;))]. Hence, 
the resonance frequencies of the incident photon are given by 
the roots of, in general, nonlinear equations u> = Re[wi(w)], 
co = Re [ 0 ) 2 ( 0 ;)]. These equations imply that the resonance 
energies Rep5i(w)], Rep^w)] (and their widths ImpA ( 0 ;)], 
Im[w 2 (o;)]) depend on the frequency of incident photon, 
which comes in (l53l > via the the wave vector k = iojv g . This is 
a general feature of non Markovian behavior when the photon 
mediated interaction between qubits is not instantaneous and 
the retardation effects have to be included. In our method the 
retardation effects are automatically included since the quan¬ 
tity kd explicitly enters the expressions for the transmission 
and reflection factors. Markovian case corresponds to long 
wavelength limit, kd « 1 when the propagation time of 
photons between the qubits can be neglected and, hence, the 
qubits interact instantaneously. 

It is important that the resonance widths are directly related 
to the lifetime of qubit superposition state ( 1671 ) since they de¬ 
fine the poles of denominator D^^uj) in the low half of the 
complex energy plane. 

Below we discuss the experimental detection of resonance 
frequencies and their widths. One of the way is to Fourier 
transform the data collected in photon-photon correlation 
measurements (26]]. However, the experiments of this type 
demand serious attention to optimizing both the measuring 
system and experimental conditions. Here we suggest to ex¬ 
tract resonance parameters from directly measured transmis¬ 
sion data. As an example we consider here two identical 
qubits. 

The resonance structure of transmission (l62l > and reflection 
(l63l > is masked by the frequency dependence of their numera¬ 
tors. This obstacle can be overcome by appropriate processing 
of the output transmission data. In order the resonance peaks 
to reveal themselves we propose to divide the transmission 
(l62l > by the factor [(w — fl)/fl] 2 . Thus, we analyze the spec¬ 
tral function S(ui), which contains pure resonance structure: 

n 2 

S ^ {uj-n + iT) 2 +T 2 e 2i % kod (70) 

The plot of S('jj) which exhibits two peaks corresponding 
to solutions of two equations (see (l55l >) 

u> = Q, ± Tsin (71) 


is shown together with transmission at Fig]4] for kod = 

7r/2,r/n = 0 . 2 . 



FIG. 4: Color online. Frequency dependence of the transmission (left 
axis, black line) and spectral function (right axis, green line)for two 
identical qubits, kod = n/2, T/fl = 0.2. 

We notice that the positions of these peaks at the frequency 
axis do not coincide with the points of the full transmission. 
The latter points are located exactly where the reflection is 
zero. This is well illustrated in Fig0 and Fig]6] where the 
transmission pattern and corresponding resonance spectrum 
are shown for kod = 5.5n, T/fl = 0.2. 



FIG. 5: Color online. Transmission pattern for two identical qubits. 

kod = 5.5 tt, T /Q = 0.2. 

The transmission pattern exhibits 12 points of the 
full transmission within the range of the frequency 
axis while there are only six resonances at the fre¬ 
quencies which are given by the roots of equation 
CLO): uj/fl = 0.805,0.866,0.929,1.070,1,133,1.194 

with the corresponding widths T/fl = 

—0.155, -0.349, -0.013, -0.013, -0.349, -0.155 Only 
four resonances which are sufficiently close to real axis 
are visible in Fig. [6] It is worth noting that there are two 
resonances with the widths being much smaller the width of 
individual qubit. Hence, for kod = 5.57T these two resonances 
give the smallest decay rates for two-qubit superposition state 
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Hence, if the photon is in resonance with the first qubit, 
the second qubit remains unexcited. If the photon is in reso¬ 
nance with the second qubit, the first qubit is unexcited only 

if il^d/Vg = 7 r. 


C. Three qubits in a waveguide 

1. Spectral properties of effective Hamiltonian 


FIG. 6: Color online. Frequency dependence of the spectral function 
for two identical qubits, kod = 5.5n, Y / fi = 0.2. The y-axis is in 
log scale 


<[67]!. 

We may conclude that as the distance between qubits, d is 
increased at fixed F, a number of the roots of ( [7T1) . which for 
small T’s lie in the narrow range fi ± T, is also increased, 
being approximately equal to kod/tr, while the widths of cor¬ 
responding peaks are decreased. The latter effect is a direct 
manifestation of non Markovian behavior. 


5. Photon wave function for two qubits in a waveguide 

Photon wave function for two qubits is calculated from (1411 
with x\ = —d/2, X 2 = +d/2 and the matrix R m ,n from d58l) . 
Outside the qubit array x > d/2, x < —d/2 the wavefunction 
is given by <f44t with t and r from (I6f)b and (f6~il . Below we 
write the photon wavefunction for two qubits in the interme¬ 
diate region {—d/2 < x < +d/2 ): 


* 2 (ar) = (t l ^ (e ikx {uj - H 2 + iV 2 ) - iT 2 e ikd e~ lkx ) 
D 2 {uj) 

(72) 

It can easily be verified that the wavefunctions (l44l > and ( [721 
are continuous at the points x = ±d/2. At resonance with the 
first qubit (w = Hi) photon is reflected from the first qubit 
and does not penetrate in the inter qubit region x > —d/2. 
However, at resonance with the second qubit (w = H 2 ) the 
wave function ^f 2 {x) ^ 0 at inter qubit region — d/2 < x < 
d/2, but \D 2 (d/2) = 0 as it follows from continuity condition. 

From (l43l we calculate the probability amplitude for the 
first or second qubit to be excited. 


<1|*Q> 


g—ikd/2 

hD 2 (uj) 


Ai (w 


Yl 2 + iT 2 )-i\ 2 ^/TJ r 2 e 2ikd 
(73) 


<2|*q> = 


gikd/2 

hD 2 {ui) 


A 2 {uj 


fit + iFi) — iAiy/TiT^ 

(74) 


Here we consider three noninteracting qubits in a waveg¬ 
uide . The qubits are positioned at the points x\ = —d, 
x 2 = +d and x 3 = 0, respectively, with a distance d be¬ 
tween adjacent qubits. The Q- subspace is formed by three 
state vectors |1) = \e 1 , g 2l g 3 , 0), |2) = \gi,e 2 ,g 3 ,0) and 
|3) = 151 , 52 , 63 , 0 ). The states |1) and |2) correspond to 
qubits located at the points x = ±d, respectively. The state 
13) is for the qubit placed at the point x = 0. The P-subspace 
is formed by the vectors | k) = I51, 52 , 53 , k). According to 
(l39t the matrix of effective Hamiltonian is as follows: 

/ £i-ihV 1 -ihVT\T^e 2ikd -ihVT/T^e ikd 

H eff= [ -ihs/T\T^e 2ikd e 2 -ihV 2 -ihy/T^e lkd 

\ -ihy/t/rie ikd -ihy/T 2 T/ i e lkd e 3 -ihV 3 

(75) 

where e, and Fj, (i = 1, 2, 3) are defined in (f34l and (l40l . 
respectively. 

The roots of this Hamiltonian in the complex frequency 
plane are defined by the equation 

det i (Hi + H 2 + fl 3 ) - R e g/hj = 0 
that can be expressed as: 


(cu — Hi + iTi) (u> — H 2 + *r 2 ) (uj — H 3 + ir 3 ) 

+ (w - Hi + *ri) r 2 r 3 e 2ifed + (w - n 2 + *r 2 ) rir 3 e 2ifcd 
+ (w - fi 3 - »r 3 ) rir 2 e 4ifed = o 

(76) 

We note that in general the energies and the widths of reso¬ 
nances depend on the physical frequency to (k = oj/v g ). 

For identical qubits (Hi = fi 2 = H 3 = H, Ti = T 2 = 
r 3 = T) we obtain from d76l > 

(u> - H + *F ) 3 + 2 (cD - H + zT) r 2 e 2ifcd , 

+ (Q - H - ?r) r 2 e 4ifcd = 0 1 " 1 

In the long wavelength limit we find from ( 1771 ) 

(w - H ) 3 + 3ir(w - H ) 2 = 0 

which gives two resonances with null width and one reso¬ 
nance which absorbs the widths of all three qubits: 


From the definition of T (l40l> we may rewrite (1741 ) as: 
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uj 1,2 = 12 ; 0J3 = 12 — 3 *r 

We obtain the same result if the running frequency cu in (1771 1 
corresponds to kd = mr, (n = 1 , 2 ....). 

The kd- dependence of real and imaginary part of three 
complex roots of equation (T771 >. u; = Rew — iT, is shown in 
Fig|7] where x = (Rew — 12) /T, y = (T — T) /T. In 3D space 
(x,y 7 z = kd) three lines of the roots are wound with a vari¬ 
able step onto a cylindrical surface which has two radii, T and 
3T/2. So that, in the projection to x. y plane these roots form 
two circles as shown in Fig| 8 ] Every point on these circles is 
merged from black and red points of Fig. [7] which belong to 
the same root. 



x 



FIG. 7: Color online. The kd- dependence of real (a) and imaginary 
(b) part of three complex roots of equation ( 1771 . uj = Red; — iT, 
where x = (Red; — fi)/T, y = (T — T)/T. The line numbers 
correspond to real and imaginary parts of three roots of Eq|77l 


FIG. 8: The projection of the three roots of Eq. [77] to x, y plane, 
where x = (Red; — f2)/T, y = (T — Tj/T. Every point on this 
graph are merged from black and red points of Fig. [7] which belong 
to the same root. 

where 

D 3 (uj) = (uj — 12i + iTi) (uj — 12 2 + 2 ) (UJ — D 3 + ^ 3 ) 

+ (w - fii + iro r 2 r 3 e 2lfed + ( u - n 2 + *r 2 ) r^^ 21 ^ 

+ (w - 12 3 - *r 3 ) Y l T 2 e iikd 

( 80 ) 

G(u>) =r 1 (u-Sia + ir 2 )(w - n 3 + iY 3 )e~ 2lkd 

+r 2 (w — 12 1 + (F-i) (uj — Q 3 + iY 3 )e 2ikd 

+r 3 (w — 12i)(u; — 12 2 ) + ir 2 r 3 (u; — fli) (81) 

+*rir 3 (ci; — 12 2 ) — 2*rir3(u; — 12 2 + *r 2 ) 

-2irir 2 (w - n 3 )e 2ikd - 2iY 3 t 2 (uj - i2i) e 2iM 

For identical qubits we obtain 


2. Transmission and reflection spectra for three qubits in a 
waveguide 


Transmission and reflection factors are calculated from (1451 
and (l46l) . where rn. n = 1,2,3 and R mn is given in Appendix. 
In the frequency picture with E = Hu; — h(Yl 3 + 12 2 + ^ 3 )/2 
we obtain for three qubits t and r the following expressions: 


(ca — 12i)(u; — 12 2 )(u; — 123) 
D 3 (uj) 


(78) 


(uj — 12 ) 

Df M 


.GidW) 

^ l Df(uj) 

where D ld (uj) is calculated in the section 2 of Appendix. 

Dl d (uj) = (uj-{ 2 + iT) 3 + 2 (w - 12 + *r) r 2 e 2ifcd 
+ (uj-n-iY) r 2 e' likd 


(83) 


(84) 


r 3 = -i 


D 3 (uj) 


(79) 


G id (uj) = 2r[(w - 12 ) 2 - r 2 ] cos 2kd + r(w - 12 ) 2 + 2r 3 
+4T 2 (w — 12) sin 2 kd 


(85) 
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In the long wavelength limit (kd << 1) we obtain from (182} 
and ( 183} 


*3 


ui — Q 

u! — ci + 3ir 


( 86 ) 


3*r 


r 3 = - 


to — ci + 3*r 


(87) 


We note that the expressions ([ 86 }. (187} are valid for a 
broad range of frequencies satisfying the condition kd « 1 . 
However, as k = uj/v g these expressions are also valid for 
kd = mr,n = 1 , 2 ,.. but only at the fixed frequencies 

uj n = nirvg/d. 

As in the two qubits case d65} . (166} the expressions ([ 86 } . 
(187} are similar to the ones for one qubit ( |49} . ([50} but with 
the width that is three times greater. 

Below we show several plots of transmission and reflection 
amplitudes for different values of kod , where fco = 

The plots are calculated for three identical qubits from (182} 
and (183} . 



FIG. 9: Color online. The dependence of transmission (black) and 
reflection (red) amplitudes on the frequency of incident photon, u>/Cl 
for different values of kgd for three identical qubits. 


3. Photon mediated entanglement for three qubits 

Analogous to two qubit case the function of the qubit sys¬ 
tem T'q (l 20 t can be written as a linear superposition of the 
three two-qubit states T'q = a |l)+6 |2)+c |3), where, in gen¬ 
eral, a, b, and c depend on the physical frequency ui. For three 
identical qubits we obtain a general expression which de¬ 
scribes the frequency dependent entanglement of three three- 


qubit states: 

^ = WT (K w - n + *r fe- ikd + TV- 
u 3 (w) 

- iT{ui - Cl + iT)e lkd - iT(uj -Cl- *T)e 3ifcd ]|l> 

(w - Clfe ikd |2) + (w - fi)[w - n + *T( 1 - e 2ifcd )]|3)) 

(88) 

In the long wavelength limit kd « 1 the maximally entan¬ 
gled superradiant state which corresponds to a coherent sym¬ 
metric superposition of three three-qubit states is formed: 

i*Q)s = -— (|1) + |2) + |3)) (89) 

w — si + 3*1 

The transmission and reflection in this case are given by the 
expressions ([86} and ( 187} . The resonance line of superradiant 
state is directly given as the line of reflection factor d87} . 

For arbitrary values of kd maximally entangled states are 
formed only for fixed values of the frequency lo. For example, 
if kd = utd/vg = n7r (n = 1,2,...) we obtain from (188} the 
expression 

<lq = ai„ - O + 3ir I(- 1 1”(I 1 > + l 2 » + l 3 d » 

where u> n = mrv g /d. 

For on resonant excitation (u> = Cl) and k^d mr, 
where fco = Cl/v g we get from d88} unentangled state s|/q = 
i (A/r) e lkod \\). In this case we observe a full reflection with 
only the first qubit being excited. 


4. Resonances in three- qubit system 

As in the two qubit case, the resonance frequencies in d78} 
and ( 179} are given by three equations u = Re[wi(u;)], ui = 
Re[w 2 (w)]’ w = Re[w 3 (w)], where wi, w 2 and w 3 are the roots 
of equation d76} . 

The denominator ( 180} can then be written as D%(u>) = 
[uj — wi(w)] [w — w 2 (w))] [w — 0 ) 3 ( 0 ;))]. Hence, the reso¬ 
nance frequencies of the incident photon are given by the 
roots of, in general, nonlinear equations uj = Re[w 1 (cu)], 
uj = Re[w 2 (w)], uj = Re[w 3 (w)]. 

Below we consider three identical qubits. Similar to the 
case of two qubits we define a spectral function by dividing 
the transmission ( 182 } by the factor [(w — fl)/H] 3 : 


S(uj) = 


Cl 3 

Df(u) 


(91) 


where D ld (uj) is given in d84t . 

Below we show the plots of transmission and spectral 
function S(u>) which exhibits peaks corresponding to the 
roots of the equations uj = Re[wi(w)](i = 1,2,3), where 
wi.wa and uj :i are the roots of the equation ( 177} . The res¬ 
onance spectrum for k$d = 7t/2, F/ST = 0.2 is shown at 
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FigOa For this case there are three resonances at the fre¬ 
quencies oj/kl = 0.8,1.0,1.2 with corresponding half widths 
r/f2 = —0.046,-0.40,-0.046. However, only two reso¬ 
nances which are closest to the frequency axis (with the width 
2r /kl = 0.08), are visible in FigJTOl 



FIG. 10: Color online. Frequency dependence of the transmission 
(left axis, blue line) and spectral function (right axis, green linelfor 
three identical qubits, kod = 7t/2, V/Q. = 0.2. 

With the increase of the inter qubit distance d, the number 
of resonances are also increased, being in the vicinity of kl 
for small F’s. The corresponding transmission pattern and the 
resonance spectrum are shown in Fiql 1 1 land Fig[l2]for kod = 
5.57r,r/H = 0.2. 



FIG. 11: Color online. Transmission pattern for three identical 
qubits, kod = 5.5n, T /kl = 0.2. 

In this case in the range uj/kl = 0.6 : 1.4 there are 13 
resonance frequencies, ten of which are seen in Fig[l2] The 
two highest peaks have the half width F/fl = 3.5 x 10“ 3 and 
r/fl = 1.63 x 10~ 2 , respectively. We note that the widths of 
some resonances, which define the decay rates of the superpo¬ 
sition state (l88l >. are much smaller than the width of individual 
qubit. 

We note that, the points of the full transition which corre¬ 
spond to the zeros of the reflection, do not coincide with the 
frequencies of the resonances. While for some frequencies 



FIG. 12: Color online. Frequency dependence of the spectral func¬ 
tion for three identical qubits, kod = 5.5n, V/kl = 0.2. The y-axis 
is in log scale. 


the corresponding values can be close to each other, in gen¬ 
eral, many reflection zeros are out of the range of resonance 
frequencies. 


5. Photon wave function for three identical qubits in a waveguide 

Photon wave function for three identical qubits is calculated 
from SB with the matrix R mn defined in ( 1A.9I >. 


(a;) = e ikx - ihTe ik \ x ~ d \ [e ikd R 13 + e~ ikd R 12 + f?i 3 ] 

_ inTe ik\x+d\ [ e ikd Ri2 + e -ikd Rii + fl 13 ] 

-iHTe ik 1*1 [e ikd R 13 + e~ ikd R 13 + f? 33 ] 

(92) 

Outside the qubit array x < — d and x > d the wavefunction 
( l92l ) is similar to (l44l) . where the transmission t and reflection 
r are given in d82l) and (l83l) . Inside the array we obtain 

3 (s) = ~ \e ikx (oj - n + ir) - iTe~ lkx e 2ikd ] 

K ’ D\ d (uj) 1 K ' J 

(93) 


for 0 < x < d, and 


d/ofo) = 

31 J Ww 

~iTe~ ikx [(w - 


e ikx {uj-n + iT) 2 +T 2 e 2ikd 
0 + ?T) + (w-0-?T)e 2ifcd ]] 


(94) 


for —d < x < 0. 

Similar to the two qubit case, here at the exact resonance 
(at = kl) the photon wavefunction does not penetrate in the 
inter qubit region. 

In the conclusion to this section we write the probability for 
the particular qubit in the array to be in excited state. From 
(l43l) we obtain: 


UI*Q> = 


(w- 


—*F (w — kl — ir) e 


n + iT) 2 e~ ikd 
3lkd -iT (w - ft + 2iT) e ikd ] 


( 95 ) 
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«*«> = S ‘ 96) 

<3|*Q> = [W - + *r (1 - e* ikd )] ( 97 ) 

We see that at resonance the first qubit only is excited. This is 
consistent with the above conclusion that at resonance photon 
does not penetrate beyond the first qubit. 

All qubit arrays considered above have a general property: 
if the photon frequency is equal to the resonance frequency of 
any qubit in the chain, the transmission signal is absent. We 
attributed this property to the destructive interference between 
the input wave and the wave which resulted from the virtual 
transitions between the qubits and the photon field in the res¬ 
onator. However, it is not clear to what extent this property can 
be attributed to uniform distribution of the qubit in the chain 
with the equal distance between adjacent qubits. The simplest 
system where we can check this property is the three qubit 
chain. We made the calculation of the transmission for three 
different qubits which are positioned at the points x± = —di, 
X2 = +d-> and 23 = 0 , respectively, with unequal distance 
between adjacent qubits. It turned out that in this case the 
transmission is similar to (f78l) with the just the same numera¬ 
tor, but different denominator, which is given in the Appendix. 
Hence, we may assume that for nonuniform qubit array the 
transmission is also zero if the input photon is at resonance 
with any qubit in the chain. 


6 . The manipulation of the photon transmission in three qubit 
chain 

Most of artificial atoms, which are used as qubits, can be 
addressed individually, so that every qubit frequency, flj, in 
the chain can be tuned from external source. Here we show 
how this property can be used to manipulate the photon trans¬ 
mission through a waveguide. As an example we consider 
three non-identical qubits, which are positioned at the points 
X\ = —d, X 2 = +d and 23 = 0, respectively, with a dis¬ 
tance d between adjacent qubits. The frequencies fli and IF 
correspond to qubits at the points 21 , 22 , respectively, while 
O3 corresponds to central qubit at the point 23. In general, if 
all three frequencies are different we obtain for long wave¬ 
length (kd << 1) transmission the plot shown in Figfl3l 
where for comparison the transmission for identical qubits are 
also shown. The qubit frequencies are directly given by zeros 
of the transmission, while two narrow nearby peaks, which lie 
between qubit frequencies correspond to the full transmission. 

Suppose now that the frequencies of the left and right qubit 
are equal (Hi = Q 2 = fi) and we manipulate the frequency of 
the central qubit, O 3 . Then we obtain the picture like the one 
shown in Fig[l4] By changing the qubit frequency O 3 we may 
move the frequency of the transmission resonance between O 
and H 3 and manipulate its width. When SI 3 becomes equal to 
O the transmission signal disappear. 



FIG. 13: Color online. Transmission for three non-identical qubits 
(blue line). For comparison the transmission for identical qubits is 
also shown (black line) 



FIG. 14: Color online. Transmission for three qubits, two of which 
are identical. Dashed (red) line: Q 3 /CI 1 = H 3 /H 2 = 1.2. Dash- 
dotted (blue) line: H 3 /H 1 = D 3 /D 2 = 2. For comparison the 
transmission for identical qubits is also shown by solid (black) line. 


D. N qubits in a waveguide 


In principle the transmission and reflection for any num¬ 
ber of qubits can be found from general expressions (1451) and 
©. For a chain of N identical homogeneously distributed 
two level atoms the analytical expression for the ID transmis¬ 
sion was found in |[29]. While in general case it is not easy to 
find analytical solutions for N qubits, nevertheless, from pre¬ 
vious calculations of the transmission for one, (l49l >. two, (l60l >. 
and three (T78l > qubits, we may guess the general structure of 
the transmission for N qubits in a waveguide: 


II (<*> - ^n) 

, _ n—1 

N D n (lo) 


(98) 
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where 


D n { uj) = det 



N 




Heff/h 


(99) 


As for the spectral properties of the effective Hamiltonian, 
we can assert that the secular equation det (E — H e ff ) has N 
poles in the low half plain of the complex energy. For identical 
qubits in the long wavelength limit there are N— 1 stable states 
with uj = VI and one resonance uj = VI — iNT which absorbs 
all the widths of individual qubits: 


The results obtained in the paper are of general nature and 
can be applied to any type of qubits. The specific properties 
of the qubit are encoded in only two parameters: the qubit 
energy VI and the rate of spontaneous emissin T. For example, 
for a superconducting flux qubit fl = y/e 1 + A 2 where e is 
an external parameter which by virtue of external magnetic 
flux, 4>x controls the gap between ground and excited states 
14^1, and the quantity A is the qubit’s gap at the degeneracy 
point (e = 0). The rate of spontaneous emission F = gA/Vl 
[4311, where g is the qubit- waveguide coupling. 


fjv 


uj — VI 

uj — Vl + iNT 


iNT 


r N = - 


uj — fl + i NT 
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VI. CONCLUSION 


Appendix 


In this paper we develop a new technique for the investiga¬ 
tion of the photon transport through multiple qubit array in a 
ID waveguide. The technique is based on the projection oper¬ 
ators formalism and non Hermitian approach, which is known 
to be a successful tool in some fields of nuclear physics and 
condensed matter. We considered in detail the one photon 
transport for two and three qubits in a waveguide, and made 
some conclusions for N qubit case. We showed that the in¬ 
teraction of qubits with a photon field results in the frequency 
dependent superposition of the qubit states. We investigated 
in detail the resonance spectra for two and three qubits and 
showed that in non Markovian case (kd » 1) the resonance 
widths, which define the decay rates of the superposition state, 
can be much smaller than the decay width of individual qubit. 

We also showed that in the long wavelength limit for uni¬ 
formly distributed array of identical qubits a coherent superra¬ 
diance state is formed with the width being equal to the sum of 
the widths of spontaneous transitions of N individual qubits. 
Within the framework of our method it is not difficult to ac¬ 
count for the decay of the qubit states to the modes other than 
the waveguide continuum. It can be done by simply adding an 
imaginary term in the qubit energy level, Vl n —> Vl n — iT' n . 

The approach developed in the paper can be easily general¬ 
ized to include the exchange interaction Hj between nearest 
neighbor qubits: 


N 

Hj = hJ2J* (4^- +1 + < + V ) (102) 

2—1 


1. The calculation of integral J(x m , x n 

In (|37| > the energies E q and Ek are the energies of incident 
photon |g) (|fc)) plus the energy of N qubits in the ground 
state. Hence, E q — Ek = h(oj q — u>k). For (l37l > we, therefore, 
have 


A 




+oo 



—oo 


— LUk + is 


(A.l) 


The main contribution to this integral comes from the re¬ 
gion ujk ~ uj q - Since u>k is the even function of k , the poles of 
the integrand ( IA. 1 b in the k plane are located near the points 
k « ±q. For an arbitrary frequency uj q that is away from the 
cutoff of the dispersion, with the corresponding wave vector 
±q , we approximate u>k around +q and —q as 


uj k ~ uJ q + {k - q) 


duJk 

dk 


= uj q + (k - q)v g (A.2) 


k=+q 


Uk ~ i>jq + (fc + q) 


diJk 

dk 


= ujq- (k + g)t 


(A.3) 


k=—q 


Near the poles the denominator in (IA.1I) takes the form: 

-v g (h-q) + i£ (A.4) 


v g (k + q) +ie 


(A.5) 


For this case it is necessary to change only the matrix of ef¬ 
fective Hamiltonian d39l >: 

(m| H eff | n) — £m&m,n “1“ n—1 $m,n—l + m,n+ 1 

- ^(r m r n ) 1 / 2 e ifc i d —i (103) 


Therefore, one pole is located in the upper half of the k plane, 
k = q + is, the other pole is located in the lower half of the k 
plane, k = —q — is. For positive d rnn , when calculating the 
integral (IA.1I) we must close the path in the upper plane. For 
negative d mn the path should be closed in lower plane. Thus, 
we obtain: 
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where 


J(a 


2 m c ik\d mn | 

hVg 


(A.6) 


2. Calculation of the f? matrix for three qubits in a waveguide 


Df{u) = {u- n + ir ) 3 4 5 6 + 2r 2 (w - n + *r) e 2ifcd 
+r 2 (w - n - *r) e 4ifed 

(A. 10) 


The matrix (m, n = 1.2.3) is calculated as the in¬ 

verse of the matrix (E — f? e //)m,n, the elements of which 
can be found from (1751 1. The matrix „ is symmetric so 
that R 12 = R 2 \,Ri 3 = R 31 , f ?23 = i? 3 2 - Direct calculations 
yield for R m ,n the following result: 


D 3 (E)R 11 

= (E~ 

£2 + 

ihY 2 ) {E - 

£3 

+ iKY 3 ) 

+h 2 Y 2 Y 3 e 2lkd 





D 3 (E)R 22 

= (E~ 

£i + 

iHY 2 ) (E - 

£3 

+ iHY 3 ) 

+h 2 Y 1 Y 3 e 2lkd 





D 3 {E)R 33 

= (E- 

£2 + 

ihY 2 ) (E - 

£1 

+ ihY 3 ) 

+h 2 Y{Y 2 e iikd 





D 3 (E)R 12 

= ~(E 

- £3 

+ iHY 3 ) ihy/Y\ 

iY 2 e 2lkd 


-h 2 ^Tj^Y 3 e 2ikd 


D 3 (E)R 13 = -(E-e 2 + *nr 2 ) 

-h 2 y/vtr 3 T 2 e 3lkd 

D 3 (E)R 23 = -(f?- ei + t^Ti) *^yr^e ifcd 

-ft 2 \/T^TV"ie 3lfcd 


3. The transmission for three qubit chain with unequal 
distance between each other 


Here we consider three different qubits which are posi¬ 
tioned at the points X\ = —d\, x 2 = +d 2 and x 3 = 0, re¬ 
spectively, with di ^ d 2 . The calculations yields the result: 


t 3 


(uj — Di)(cu — Q 2 )(ut — Q 3 ) 


(A.ll) 


where D 3 (E) = det (E - H eS ) mn : 


where 


Ds(E) — (E — E\ + iKF 1 ) (E — £2 + ih ^* 2 ) (E — £3 + ihY 3 ) 

+ (E-e 1 + iKTi) H 2 T 2 T 3 e 2ikd + 

(E-e 2 + ihV 2 ) h 2 T 1 Y 3 e 2ikd + (E - e 3 - ihY 3 ) H 2 Y 1 Y 2 e 4ikd 

(A. 8 ) 

and the quantities Ei(i = 1, 2,3) are defined in (l34l >. 

At the end of this subsection we write down from (IA. 7b the 
matrix R m ,n for three identical qubits. 

hD l 3 d {u>)Ru = (u- fi + IT ) 2 + r 2 e 2ifed 

HD l 3 d (uj)R 33 = (w - fi + IT ) 2 + r 2 e 4ifed 

HD l 3 d (uj)R 12 = -iY (w - fi) e 2lfcd (A.9) 

HDf (J)R 13 = -iY (u - n + iT) e ifed - T 2 e 3ifcd 

i?22 = Rll,E,23 = Rl3 = R w _ 7791 = -R 19 _ 


F{uj) = C/if/ 2 t/ 3 + iY ± Y 2 Y 3 (e 2ikdl + e 2ikd2 - e 2ifc ( d i +d2 ) - l) 
+t/3riT 2 ( e 2ikd '+ 2ikd * - 1) + U!Y 2 Y 3 (e 2ikd2 - l) 

+[/ 2 rir 3 (e 2ifcdl - 1) + iUiU 2 Y 3 + ilhU 3 Y 2 + iU 2 U 3 Y 

(A. 12) 

U\ — u> — fiij U 2 = to — fl 2 ; U 3 = cj — Yl 3 

The equation (IA. 1 lb is similar to (f78l> with the same numer¬ 
ator but different denominator. 


[1] I. Buluta, S. Ashhab, F. Non. Natural and artificial atoms for 
quantum computation. Rep. Progr. Phys. 74. 104401 (2011). 

[2] J.Q. You, F. Nori, Atomic physics and quantum optics using 
superconducting circuits. Nature 474, 589 (2011). 

[3] I. Buluta, F. Nori, Quantum Simulators. Science 326, 108 
(2009). 

[4] Yao-Lung L. Fang and H. U. Beranger, Wavegude QED: Power 
spectra and correlations of two photons scattered off multiple 
distant qubits and a mirror. Phys. Rev. A91, 053845 (2015). 

[5] A. Blais, R-S. Huang, A. Wallraff, S. M. Girvin and R. J. 
Schoelkopf. Cavity quantum electrodynamics for supercon¬ 
ducting electrical circuits: an architecture for quantum com¬ 
putation Phys. Rev. A 69, 062320 (2004). 

[6] A. Wallraff, D. I. Schuster, A. Blais, L. Frunzio, R-S. Huang, 


J. Majer, S. Kumar, S. M. Girvin and R. J. Schoelkopf, Strong 
coupling of a single photon to a superconducting qubit using 
circuit quantum electrodynamics. Nature 431 162 (2004). 

[7] R. J. Schoelkopf and S. M. Girvin, Wiring up quantum systems. 
Nature, 451, 664 (2008). 

[8] Astafiev O, Zagoskin A M, Abdumalikov A A, Pashkin Yu A, 
Yamamoto T, Inomata K, Nakamura Y and Tsai J S., Reso¬ 
nance fluorescence of a single artificial atom. Science 327, 840 
( 2010 ). 

[9] Astafiev O V, Abdumalikov A A, Zagoskin A M, Pashkin Yu A, 
Nakamura Y and Tsai J S Ultimate on-chip quantum amplifier. 
Phys. Rev. Lett. 104, 183603 (2010). 

[10] Abdumalikov A A, Astafiev O, Zagoskin A M, Pashkin Yu A, 
Nakamura Y and Tsai J S Electromagnetically induced trans- 












16 


parency on a single artificial atom Phys. Rev. Lett. 104 ,193601 
( 2010 ). 

[11] Abdumalikov A A. Astafiev O V, Pashkin Yu A, Nakamura Y 
and Tsai J S.. Dynamics of coherent and incoherent emission 
from an artificial atom in a ID space. Phys. Rev. Lett. 107, 
043604 (2011). 

[12] Hoi I-C, Wilson C M, Johansson G, Palomaki T. Peropadre B 
and Delsing P Demonstration of a single- photon router in the 
microwave regime Phys. Rev. Lett. 107, 073601 (2011). 

[13] Hoi I-C. Palomaki T, Lindkvist J, Johansson G. Delsing P and 
Wilson C M Generation of nonclassical microwave states us¬ 
ing an artificial atom in ID open space. Phys. Rev. Lett. 108, 
263601 (2012). 

[14] I.-C. Hoi, A. F. Kockum, T. Palomaki, T. M. Stace, B. Fan. 
L. Tornberg, S. R. Sathyamoorthy, G. Johansson. P. Delsing, 
and C. M. Wilson, Giant cross Kerr effect for propagating mi¬ 
crowaves induced by an artificial atom. Phys. Rev. Lett. Ill, 
053601 (2013). 

[15] Hoi I-C, Wilson C M, Johansson G, Lindkvist J, Peropadre B. 
Palomaki T and Delsing P Microwave quantum optics with an 
artificial atom in one-dimensional open space. New J. Phys. 15, 
025011 (2013). 

[16] A. F. van Loo, A. Fedorov, K. Lalumiere, B. C. Sanders, A. 
Blais, and A. Wallraff, Photon-mediated interactions between 
distant artificial atoms. Science 342, 1494 (2013). 

[17] L. DiCarlo, M. D. Reed, L. Sun, B. R. Johnson, J. M. Chow, 
J. M. Gambetta, L. Frunzio, S. M. Girvin, M. H. Devoret, and 
R. J. Schoelkopf, Preparation and Measurement of Three-Qubit 
Entanglement in a Superconducting Circuit. Nature 467, 574 
( 2010 ). 

[18] N. Roch, M. E. Schwartz, F. Motzoi, C. Macklin, R. Vijay, A. 
W. Eddins, A. N. Korotkov, K. B. Whaley, M. Sarovar, and 
1. Siddiqi, Observation of measurement-induced entanglement 
and quantum trajectories of remote superconducting qubits. 
Phys. Rev. Lett. 112, 170501 (2014). 

[19] M. Jerger, S. Poletto, P. Macha, U. Hubner, A. Lukashenko, E. 
Il’ichev and A. V. Ustinov Readout of a qubit array via a single 
transmission line. Eur. Phys. Lett. 96, 40012 (2011). 

[20] M. Jerger, S. Poletto, P. Macha, U. Hubner, E. ITichev, and 
A. V. Ustinov, Frequency Division Multiplexing Readout and 
Simultaneous Manipulation of an Array of Flux Qubits. Appl. 
Phys. Lett. 101, 042604 (2012). 

[21] P. Macha, G. Oelsner, J.-M. Reiner, M. Marthaler, S. Andre, G. 
Schon, U. Huebner, H.-G. Meyer, E. Il’ichev, A. V. Ustinov, Im¬ 
plementation of a Quantum Metamaterial using superconduct¬ 
ing qubits. Nat. Commun. 5, 5146 (2014). 

[22] J.-T. Shen and S. Fan. Coherent Single Photon Transport in 
a One-Dimensional Waveguide Coupled with Superconducting 
Quantum Bits. Phys. Rev. Lett. 95, 213001 (2005). 

[23] J.-T. Shen and S. Fan, Coherent photon transport from sponta¬ 
neous emission in one-dimensional waveguides. Optics Letters 
30, 2001 (2005). 

[24] Yao-Lung L Fang, Huaixiu Zheng and Harold U Baranger, One¬ 
dimensional waveguide coupled to multiple qubits: photon- 
photon correlations. EPJ Quantum Technology 1, 3, (2014). 

[25] J.-T. Shen and S. Fan, Theory of single photon transport in a 
single-mode waveguide. I. Coupling to a cavity containing a 


two-level atom. Phys. Rev. A79, 023837 (2009). 

[26] H. Zheng and H. U. Baranger, Persistent Quantum Beats and 
Long-Distance Entanglement from Waveguide-Mediated Inter¬ 
actions. Phys. Rev. Lett. 110, 113601 (2013). 

[27] K. Lalumiere, B. C. Sanders, A. F. van Loo, A. Fedorov, A. 
Wallraff, and A. Blais, Input-output theory for waveguide QED 
with an ensemble of inhomogeneous atoms. Phys. Rev. A88, 
043806 (2013). 

[28] S. Fan, S. E. Kocabas, and J.-T. Shen, Input-output formal¬ 
ism for few-photon transport in one-dimensional nanophotonic 
waveguides coupled to a qubit. Phys. Rev. A 82, 063821 (2010). 

[29] T. S. Tsoi and C. K. Law, Quantum interference effects of a 
single photon interacting with an atomic chain. Phys. Rev. A78, 
063832 (2008). 

[30] Y. Xu, Y. Li, R. K. Lee, and A. Yariv, Scattering theory analysis 
of wavegude-resonator coupling. Phys. Rev. E62, 7389 (2000). 

[31] G.Diaz-Camacho, D. Porras, and J. J. Garcia-Ripoll, Photon- 
mediated qubit interactions in ID discrete and continous mod¬ 
els. Phys. Rev. A 91, 063828 (2015). 

[32] H. Feshbach, Unified Theory of Nuclear Reactions. Ann. Phys. 
(N.Y.) 5, 357 (1958). 

[33] V. V. Sokolov and V. G. Zelevinsky, Collective Dynamics of 
Unstable Quantum States. Ann. Phys. (NY) 216, 323 (1992). 

[34] S. Sorathia, F.M. Izrailev, G.L. Celardo, V.G. Zelevinsky, and 
G.P Berman, Internal chaos in an open quantum system: From 
Ericson to conductance fluctuations. Eur. Phys. Lett. 88, 27003 
(2009). 

[35] G. L. Celardo, A. M. Smith, S. Sorathia, V. G. Zelevinsky, R. 
A. Sen’kov, and L. Kaplan, Transport through nanostructures 
with asymmetric coupling to the leads. Phys. Rev. B 82, 165437 
( 2010 ). 

[36] Ya. S. Greenberg, N. Merrigan, A. Tayebi, and V. Zelevinsky, 
Quantum signal transmission through a single-qubit chain. Eur. 
Phys. J. 86, 368 (2013). 

[37] N. Auerbach and V. Zelevinsky, Super-radiant dynamics, door¬ 
ways, and resonances in nuclei and other open mesoscopic sys¬ 
tems. Rep. Progr. Phys. 74, 106301 (2011). 

[38] E. Ilichev, S. N. Shevchenko, S. H. W. van der Ploeg, M. Graj- 
car, E. A. Temchenko, A. N. Omelyanchouk, and H.-G. Meyer, 
Multiphoton excitations and inverse population in a system of 
two flux qubits. Phys. Rev. B 81. 012506 (2010). 

[39] E. A. Temchenko, S. N. Shevchenko, and A. N. Omelyanchouk, 
Dissipative dynamics of a two-qubit system: Four-level lasing. 
Phys. Rev. B 83, 144507 (2011). 

[40] S. Zippilli, M. Grajcar, E. Ilichev, and F. Illuminati, Simulating 
long-distance entanglement in quantum spin chains by super¬ 
conducting flux qubits. Phys. Rev. A 91, 022315 (2015). 

[41] I. Rotter, Dynamics of quantum systems. Phys. Rev E 64, 
036213 (2001). 

[42] Caspar H. van der Wal, A. C. J. ter Haar, F. K. Wilhelm, R. N. 
Schouten, C. J. P. M. Harmans, T. P. Orlando, Seth Lloyd, and 
J. E. Mooij, Quantum Superposition of Macroscopic Persistent- 
Current States. Science 290, 773 (2000). 

[43] A. N. Omelyanchouk, S. N. Shevchenko, Ya. S. Greenberg, O. 
Astafiev, and E. Il’ichev, Quantum behavior of a flux qubit cou¬ 
pled to a resonator. Low Temp. Phys. 36, 893 (2010). 



